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Abstract 



O 

<N : The non-equatorial spherical null geodesies of rotating Kerr black holes are studied analyti- 



cally. Unlike the extensively studied equatorial circular orbits whose radii are known analytically, 



O: 

■ no closed-form formula exists in the literature for the radii of generic (non-equatorial) spherical 

geodesies. We provide here an approximate formula for the radii rph(a/M; cos i) of these spherical 

null geodesies, where a/M is the dimensionless angular- momentum of the black hole and cosz is 

an effective inclination angle (with respect to the black-hole equatorial plane) of the orbit. It is 

well-known that the equatorial circular geodesies of the Kerr spacetime (the prograde and the 

^sO . retrograde orbits with cosi = ±1) are characterized by a monotonic dependence of their radii 
OO 

'nJ" ■ rph(a/M; cos i = ±1) on the dimensionless spin-parameter a/M of the black hole. We use here our 
CN ■ 

■ novel analytical formula to reveal that this well-known property of the equatorial circular geodesies 
, is actually not a generic property of the Kerr spacetime. In particular, we find that counter-rotating 

spherical null orbits in the range (3\/3 — \/59)/4 ^ cosi < are characterized by a non- monotonic 

^ ■ dependence of rph(a/M;cosi = const) on the dimensionless rotation-parameter a/M of the black 

5t ■ hole. Furthermore, it is shown that spherical photon orbits of rapidly-rotating black holes are 



characterized by a critical inclination angle, cosi = y 4/7, above which the coordinate radii of the 
orbits approach the black- hole radius in the extremal limit. We prove that this critical inclination 
angle signals a transition in the physical properties of the spherical null geodesies: in particular, it 
separates orbits which are characterized by finite proper distances to the black-hole horizon from 
orbits which are characterized by infinite proper distances to the horizon. 
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I. INTRODUCTION 



The characteristic geodesies of the Kerr black-hole spacetime have been extensively stud- 
ied since the pioneering work of Carter see also {2 7| and references therein. Of particular 
importance are spherical null orbits — orbits with constant coordinate radii on which mass- 
less particles can orbit the black hole. The spherical null geodesic ( "photonsphere" ) of a 
black hole provides valuable information on the structure and geometry of the black-hole 
spacetime. 

The spherical null orbits are especially interesting from both an astrophysical and theo- 
retical points of view. For example, the optical appearance to external observers of a star 
undergoing gravitational collapse is related to the physical properties of the photonsphere 



fiai 



|9| . This surface also determines how the night sky would appear to an observer near a 
black hole or a very compact star [lOi]. In addition, spherical null orbits are closely related 
to the characteristic scattering resonances of black holes in the eikonal limit (the geometric- 

n nn 

optics approximation) [6|, |11|-|16|. According to the geometric-optics approximation, these 
characteristic black-hole quasinormal resonances correspond to massless particles trapped 



at the unstable null orbit and slowly leaking out 
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11 



16|. 



Furthermore, it was recently proved that for hairy black-hole configurations, the black- 



hole photo,.phe.e provides a genenc lowe. bound on the effective length of the hait 

In addition, it was recently proved \19\ that circular null geodesies provide the fastest way 
to circle a black hole as measured by asymptotic observers. 



It is worth emphasizing that in most cases of physical interest [20|, |2l| the physical 
properties of the spherical null geodesies (and in particular, the radii and the energy-to- 
angular-momentum ratio) must be computed numerically. The only known exceptions are 
the prograde and retrograde circular orbits in the equatorial plane of the rotating Kerr black 
hole, in which case closed analytical formulae for the physical properties of the orbits have 
been provided in see Eq. ( 12T]) below. The parameters of the null polar orbit are also 



known in a closed form [Sl], see Eq. ( l22l) below. 

To the best of our knowledge, no closed-form formula exists in the literature for the radii 
of generic (non-equatorial) spherical null geodesies of rotating Kerr black holes. One of the 
main goals of the present paper is to derive effective analytical formulae which describe the 
physical properties of such non-equatorial spherical geodesies in the Kerr spacetime. 
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It is well-known that the equatorial (prograde and retrograde) circular orbits of the 
Kerr spacetime are characterized by a monotonic dependence of their radii rph(a) on the 
spin-parameter a of the black hole [2']. In the present study we shall use our new analytical 
formula [see Eq. ([21]) below] to show that this well-known property of the equatorial circular 
geodesies is actually not a generic property of the Kerr geometry. In particular, we shall 
show that there is a significant fraction of counter-rotating spherical orbits for which the 
function rph(a) exhibits a non- monotonic behavior. 



Recently, Yang et. al. 16| have studied the null spherical orbits of the Kerr spacetime. 



One of the most remarkable conclusions of 



16| is that near-extremal Kerr black holes are 



characterized by a significant fraction of spherical null geodesies whose radii approach the 
)lack-hole radius in the near-extremal limit. In particular, it was observed numerically in 
1^ that such near-horizon orbits exist for near-extremal black holes in the finite interval 

sin = 0.731 < sin0 < 1 , (1) 

where 9 is the inclination angle of the orbit with respect to the polar axis. Below we shall 
provide a fully analytical explanation for this phenomena. Furthermore, we shall obtain an 
analytical expression for the exact value of the critical polar angle sin 6'c above which the 
near-horizon spherical null geodesies appear. 

The rest of the paper is devoted to the investigation of the physical properties of non- 
equatorial spherical null geodesies in the rotating Kerr spacetime. In Sec. II we describe 
the dynamical equations which determine the null geodesies of the Kerr spacetime. In 
Sec. Ill we obtain the characteristic equation which determines the radii of the spherical 
null geodesies. In Sec. IV we solve the characteristic equation using a series expansion of 
Tph (a/M; cos i) in powers of the dimensionless rotation-parameter a/M of the black hole. 
We then use this analytic approach to reveal the regime in which the function rph(a/M) 
exhibits a non-monotonic behavior. In Sec. V we analyze the spherical null geodesies of 
rapidly-rotating black holes and discuss the near-horizon orbits of these near-extremal black 
holes. We conclude in Sec. VI with a summary of the main results. 
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II. DESCRIPTION OF THE SYSTEM 



We shall analyze the spherical null geodesies which characterize the Kerr black-hole space- 
time. In Boyer-Lindquist coordinates the metric is given by (we use gravitational units in 



which G = c=l) [3|,l22| 



2Mr 



, 2 AMar sin^ 6^ , , , p' 2 
at ata(p + —ar 



+p'de' + (r' + a' + 



2MaV sin^ 6 
? 



sm 



,2 

a' 

2 nji2 



(2) 



where M and a are the mass and angular momentum per unit mass of the black hole, 
respectively. Here A = — 2Mr + and p = + cos^ 6. The black-hole (event and 
inner) horizons are located at the zeroes of A: 



r± = M ± (M^ - a 



2x1/2 



(3) 



Null 



motion 



geodesies in the rotating Kerr spacetime are characterized by three constants of the 
In terms of the covariant Boyer-Lindquist components of the 4- momentum, these 
conserved quantities are 

E = —pt = total energy , (4) 
Lz = P(p = component of angular momentum parallel to the symmetry axis , (5) 

and 

Q^pI + cos' e[-aYt + pV sin' 9] . (6) 
The geodesies in the black-hole spacetime are governed by the following set of equations 



p-Y = (LJ sin' e-aE) + aT/A , 
= a(Lz - aE sin' 9) + ir' + a')T / A 

dX 



(7) 
(8) 
(9) 
(10) 



where A is an afiine parameter along the null geodesies. Here 



T = E{r' + a') - L^a 



(11) 



4 



Vr = T^-A[{L,-aEf + Q], (12) 

Ve = Q- cos^ e[-a^E^ + Ll/ sin^ 6] . (13) 

The effective potentials Vr and Vg determine the orbital motions in the r and 6 directions, 
respectively. 

Equatorial orbits are characterized by Q = 2|. It is convenient to quantify the deviation 
of a generic (non-equatorial) orbit from the equatorial plane of the black hole using an 



effective inclination angle i which is defined by 23l-|26| 



cos 2 



L 



(14) 



where 



L ^ y^l4 + Q . (15) 

Note that L and i are constants of the motion. For spherical black-hole spacetimes (with 
a = 0), L is the total angular momentum of the particle The extensively studied simple 
equatorial geodesies are characterized by cos^i = 1, where cosi = +1/— 1 correspond to 
prograde/retrograde orbits, respectively. Polar orbits are characterized by cosi = 0. 



III. SPHERICAL NULL ORBITS IN THE BLACK-HOLE SPACETIME 

Spherical geodesies in the black-hole spacetime are characterized by the two conditions 

m 

Vr = and V'.=0. (16) 
Substituting (fT2l) into (fT6i) . one finds the ratio 

— a"^ sm t 

and the characteristic equation 

- 3Mr^ + 2aV^ sin^ i - 2MaV^ + a'^r sin^ i + Ma^ sin^ i 

+2Mar cos 2-^3r4 + (1 - 3 sin^ i)a'^r^ - sin^ i = , (18) 

which determines the radii rph(a/M; cosi) of the null spherical geodesies. It is convenient to 
define the dimensionless radial coordinate f and the dimensionless spin-parameter d of the 
black hole: 

f = r/M and d = a/M, (19) 
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in terms of which the characteristic equation (fT8|) becomes 



-5 0-4 I 0-2-3 • 2 • 0-2-2 , -4- • 2 • , -4 • 2 • 

r — 6r + 2a r sm t — 2a r + a r sm t + a sm z 



+2afcosiy3f4+(l-3sin^i)a2f2-a4sin2^ = . (20) 

The exact (analytical) solution of the characteristic equation ( l20l) is only known for the 
following two limiting simple cases: 

• The equatorial null geodesies. As mentioned above, these orbits are characterized by 
cosi = ±1, where the upper sign corresponds to the prograde circular geodesic while 
the lower sign corresponds to the retrograde circular geodesic. In this simple case one 
finds 



(a) = Tnh ( a 



,2 



■phl"J = 'phi"; cosi = ±1) = 2{l + cos[- cos ^(=Fa)]}- (21) 

The polar null geodesic. As rnentioned above, this orbit is characterized by cosi = 0. 
In this simple case one finds 



r°h(a) = rph(a; cosi = 0) = 1 + 2^1 - ^a^ cos cos ^ ^^^^3 }• (22) 

To the best of our knowledge, no closed-form formula exists in the literature for the radii 
of generic (cosz 7^ 0, ±1) spherical null geodesies of the rotating Kerr spacetime. Below we 
shall derive approximate analytical formulae for the radii of these important orbits. 

IV. EXPANSION IN POWERS OF a/M 

We shall first consider slowly-rotating black holes. Substituting the power series 

00 

fph(a;cos2) = ^^Cfc(cosi) ■ a^ (23) 

A:=0 

into the characteristic equation ( 120|) . one finds the analytical formula 

9 1 

fph(a; cos i) = 3 ^ cos i ■ a + - (— 4 + 2 cos i) ■ a 

V 3 9 

+ ^(cos=^z-2cosz) -03 + 0(61^) (24) 
81 

for the radii of the null spherical geodesies. 



In order to test the accuracy of the approximate formula i^^, we present in Table [T] 
the dimensionless ratio r'^y^^"^ /r'^'^^ between the approximate radii as obtained from the 
analytical formula (I24p and the exact radii which are obtained by numerically solving the 
characteristic equation (120]) 27l |. Although the expansion of fph(a;cosi) in powers of a is 
formally valid for slowly-rotating black holes, one learns from Table [T] that the agreement 
between the numerical data and the analytical formula is actually very good in almost 
the entire range of the dimensionless spin-parameter a. As is evident from Table [H the 
analytical approximation (12^ fails only for co-rotating orbits of rapidly-rotating black holes 
with a ^ 0.9. Below we shall obtain an alternative (and very accurate) analytical formula 
for the radii of spherical null geodesies of these near-extremal black holes. 



a/M 


cos i = —1 


cosz = — ^ 


cos i = 


cos i = +^ 


cos i = +1 





1 


1 


1 


1 


1 


0.2 


1.0000287 


1.0000364 


1.0000448 


1.0000509 


1.0000444 


0.4 


1.0003858 


1.0005469 


1.0007650 


1.0010268 


1.0009461 


0.6 


1.0016708 


1.0025501 


1.0043633 


1.0070245 


1.0069312 


0.8 


1.0045868 


1.0076324 


1.0168252 


1.0349897 


1.0376514 


0.9 


1.0068806 


1.0120146 


1.0312513 


1.0807893 


1.0930592 


1 


1.0098488 


1.0181206 


1.0585458 


1.2674558 


1.5161605 



TABLE I: Spherical null geodesies of rotating Kerr black holes. We display the dimensionless ratio 
^approx^_^exact ]-,gt^ggn the approximate radii as obtained from the analytical formula (|24p and the 



exact radii which are obtained by numerically solving the characteristic equation (|20p [27j . We 
display results for orbits with cosi = — 1 (the retrograde equatorial orbit), cos? = — cost = 
(the polar orbit), cosi = and cosi = +1 (the prograde equatorial orbit), and for various 
values of the dimensionless spin-parameter a = a/M of the black hole. The agreement between 
the numerical data and the analytical formula (|24p is very good in almost the entire range of the 
dimensionless rotation-parameter a. There is, however, a large discrepancy between the numerical 
results and the analytical formula ()24p for co-rotating orbits of rapidly-rotating black holes with 



a>0.9. 



It is well-known that the equatorial circular geodesies of the Kerr spacetime are character- 
ized by a monotonic dependence of fpii(a;cosi = ±1) on the dimensionless spin-parameter 



a of the black hole. Specifically, the equatorial prograde/retrograde null circular orbit (with 
cosi = +1/— 1) is a monotonic decreasing/increasing function of the dimensionless spin pa- 
rameter a, see Eq. f l2ip and also Figure 1 of {2!. We shall now use our novel formula (I24p to 
reveal that this is actually not a generic property of the Kerr spacetime. 

Inspection of Eq. fl24p reveals that, for co- rotating orbits with a given positive value of 
the conserved inclination angle (cosi > 0) the expansion coefficients are all negative (except 
from Co = 3). This implies that the radii fph(a; cosz = const) of such co-rotating orbits are 
all monotonic decreasing functions of the dimensionless spin-parameter a of the black hole 



28| . However, one finds from a second-order expansion of rph(a) [see Eq. ^^^] that counter- 



rotating orbits in the range (3-\/3 — v^59)/4 ^ cosi < are characterized by a non-monotonic 
dependence of rph(a;cosi = const) on the dimensionless spin-parameter of the black hole. 
In particular, we find that the function fph(a; cosz = const) has a global maximum at 

, ., 3-\/3cos2 

flmax cos? ~ ■ , ^ — (25) 

—4 + 2 cos^ I 

with 

_ , , 3 cos^ i . , 

?^max(cOS«) ~ 3 + ^ . (26) 

4 — 2 cos"^ I 

Equation (!26|) is very accurate even for rapidly-rotating black holes. For example, for cosi = 
—0.55 one finds from fl26l) fmax — 3.27. This value agrees to better than 1% with the 
numerically computed value of ~ 3.3 29 1. 



V. RAPIDLY-ROTATING (NEAR-EXTREMAL) BLACK HOLES 

We have seen that the series expansion ^2M provides a good approximation for the radii 
of the spherical null geodesies in almost the entire range of the dimensionless spin-parameter 
a of the black hole (see Table H]). However, this simple formula fails for co-rotating orbits 
of rapidly-rotating black holes with a ^ 0.9. We shall now derive an alternative analytical 
expression for the radii of the spherical null geodesies of these near-extremal black holes. In 
particular, we shall show that in the near-extremal limit the function fph(a; cosi) undergoes 
a "phase transition" at the critical inclination angle cosi = a/4/7 . 

For rapidly-rotating (near-extremal) black holes we define 

f± = 1 ± e and fph = 1 + 5ph , (27) 
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where e, 5ph ^ 1. From (1271) one finds 



a = 1 - eV2 + 0(e^ 



Substituting (1271) and (l28l) into the characteristic equation (120|) . one obtains 



(28) 



4 cos^ i(3 cos^ i - 4) (7 cos^ i - 4)6^ + {61 cos^ z - 220 cos^ i + 176 cos^ z - 64)5^ 

+ 16 cos^ i(4 - 3 cos2 t)e^ + 0{6\ e^) = . (29) 

The quahtative behavior of 5ph(e; cosz) in the near-extremal e — )■ limit depends on whether 
the coefficient of the 0(5^) term in Eq. ( 129|) is positive, negative, or zero. Note that this 
coefficient vanishes at 

^ = cosi - ^/I/7 = . (30) 



The solution of the characteristic equation ( l29l) is given by [30 1 



5ph(e; cosz) 



4 cos^ j 



wi34^ + 0(eVa 

2-1/3,2/3^0(^4/3) 
73/ 



for cos i - a/477 > e^/^ ; 

for - e^/^ < cos i - a/477 < 6^/3 ; 



^I|i(y477-cosi) + 0(eVe') for e^/s ^ ^4/7 _ ^ ^ ^ _ 



(31) 

From ( 13T|) one learns that the solution 5ph(e; cosi) of the characteristic equation (129|) exhibits 
a "phase transition" [from a (5ph(e — t- 0) — ^ behavior to a finite 5ph(e — )■ 0) behavior], which 
occurs in the extremal e — )■ limit at the critical inclination angle cos i = a/4/7. 



We shall now prove that the critical inclination angle, cost = a/4/7, signals a transition 
from finite to infinite proper distances of the null orbits from the black-hole horizon. The 
proper radial distance between the black-hole horizon (at = 1 + e) and the intersection 
point of a null spherical orbit (of radius fph = 1 + 5ph) with the equatorial plane of the black 
hole is given by [31I 

Ai = VS^ - e2 + In (5 + VS^ - e^) - Ine . (32) 
Substituting into ( l32l) the value of 5ph(e; cosi) from Eq. (13T1) . one obtains 

In ('2cosi+V4-3cos^ ^ _^ Qf^^^ • _ ^4/7 ^ ^2/3 

for - e^/^ < cos z - a/477 < e^/^ 
-lne + 0(lnO for e^/^ < a/477 - cos i < 1 , 



Ai{e ^ 0) = < 



\/7 cos-^ i— 4 



-ilne + 0(l) 



(33) 
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in the extremal e — )■ limit. We therefore find that the critical inclination angle, which 
occurs at cost = a/4/7 in the extremal limit, separates null orbits which have finite proper 
distances to the black-hole horizon from null orbits which have infinite proper distances to 
the horizon [32,]. 

Each spherical „ull orb.t . bounded ■„ some strip of the polar angle 0, where 

9 =71 — 9^ [16]. The two polar turning-points are determined from the requirement 
1^(6*^) = 0, see Eq. f|T3|) . We shall now evaluate the polar boundaries {9^ , 9^} of the near- 
horizon orbits discussed above (the orbits whose coordinate radii approach the black-hole 
horizon in the extremal limit and which have finite proper distances to the horizon in this 
limit). Taking cognizance of Eqs. ( IT71) and (13T]) . one obtains the dimensionless ratio 

L 



ME cosi 



4 cos"' I 



(34) 



for the null spherical orbits in the regime cosi - ^4/7 > e2/3. Substituting (El]) into (dl 
with Ve = 0, one finds that the polar turning-points (which characterize the maximal polar- 
deviation of a near-horizon orbit from the equatorial plane of the black hole) are given by 



33| 



2 „+ cos^ i — 4+ y/—15 cos"^ i + 8 cos^ i + 16 _ , 

cos 9 = ^1 . (35) 

2 cos2 i ^ ' 

For the equatorial circular orbit with cosz = 1 one finds cos^6'^ = [9'^ = |) as expected. 



For the critical orbit with cosi = a/4/7 one finds from ([35 



sin^^ = Y4-2V3 . (36) 

This analytically derived value is in excellent agreement with the numerical value of ~ 0.731 
computed in [16 1. 



VI. SUMMARY 



We have studied analytically the non-equatorial spherical null geodesies of rotating Kerr 
black holes. The main results presented in this paper are: 

(1) We have obtained the analytical expression ( 12^ for the radii fph(a; cosi) of generic (non- 
equatorial) spherical null orbits in the Kerr spacetime. This analytical formula is formally 
valid for slowly-rotating black holes. However, we have demonstrated explicitly (see Table 
H]) that the analytical formula (|24|) actually agrees very well with the exact (numerically 
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computed) values of the orbital radii in almost the entire range of the dimensionless spin- 
parameter a of the black hole. 

2) The extensively studied equatorial circular geodesies of the Kerr black hole are known 
2| to be characterized by a monotonic dependence of their radii fph(a;cosz = ±1) on the 
rotation-parameter a of the black hole. Our analysis revealed that this is actually not a 
generic property of the Kerr spacetime. In particular, it was shown that counter-rotating 
orbits in the range (3v^ — v^59)/4 ^ cosz < are characterized by a non-monotonic 
dependence of fph(a; cosi = const) on the dimensionless spin-parameter a of the black hole, 
see Eqs. ( 125|) and ( l26l) . Likewise, counter-rotating orbits in the range 3-\/3 — -\/3T ^ cosi < 
are characterized by a non-monotonic behavior of the dimensionless angular-momentum-to- 
energy ratio £ph(o; cosi = const) which is a constant of the motion. 

(3) We have obtained the analytical expression fl3T]) for the radii fph(a; cosz) of the spherical 
null geodesies of rapidly-rotating (near-extremal) Kerr black holes. This analytical formula 
reveals that the spherical photon orbits of rapidly-rotating black holes are characterized by 
a critical inclination angle, cosi = a/4/7, above which the coordinate radii of the orbits 
approach the black-hole radius in the extremal limit. We have shown that this critical 
inclination angle is the boundary between spherical null orbits which have finite proper 
distances to the black-hole horizon and spherical null orbits which have infinite proper 
distances to the horizon. Our analytical expression fl36|) for the critical (maximal) polar 



angle of such near-horizon orbits, sin e, = ^4 - 2^3, is remarkably close to the numerically 
computed 13] value sin^^^ = 0.731. 
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the two remaining (non-negligible) terms must be positive while the other non-negligible term 
must be negative. Specifically, there are three distinct regimes: (1) For e > e^/^ the 0{6^) 
term is negligible as compared to the other two terms in Eq. p9p . (2) For — e^/'^ <^ <^ e'^^'^ the 
0(^(5^) term is negligible as compared to the other two terms in Eq. (j29p . (3) For ^ <^ — g2/3 
the O(e^) term is negligible as compared to the other two terms in Eq. (j29p . 

[31] T. Jacobson, Class. Quant. Grav. 28, 187001 (2011). 

[32] Note that one finds A£ = ^ln3 for the prograde equatorial orbit with cosi = 1, in accord 
with the result of ^ . 

[33] It is worth emphasizing again that the effective inclination angle i is a constant of the motion, 
while the polar angle 9 varies along the trajectory within the interval [6^,6^], where the two 
boundaries are determined by the relation (j35p . 
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